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The approach of a parton plasma at future heavy ion colliders towards kinetic 
and chemical equilibrium is considered. A plasma with a self-consistent evolving 
parton-parton interaction strength is shown to equilibrate better and faster than 
the usual but inconsistent one with a fixed strength. We explain why as a conse- 
quence of this, a parton plasma is a unique kind of many-body system. Because 
our time evolution scheme does not require the plasma to be in either kind of 
equilibrium from the outset, out-of-equilibrium effect on particle productions can 
be revealed. We show this on photon production and discuss the implications on 
photon as a signal to detect the quark-gluon plasma. IASA 98-3, UA/NPPS-98-13 



1 Introduction 

The trophy of the game of relativistic heavy ion collisions at the future Rel- 
ativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC) is 
well known, namely the quark-gluon plasma. To reach this goal, one inter- 
rogates the particles coming out from the beginning till the very end of the 
collisions. Based on the information thus obtained, one determines whether 
a phase transition into the predicted nearly free quarks and gluons state has 
occurred. In order to perform the necessary analysis, qualitative understand- 
ing and quantitative control are very important. In this talk, we present some 
recent development in the former. We will look at the equilibration of the 
parton plasma in the following sections to determine the state of the system 
just before the phase transition. One of the reasons being that if the system 
is able to complete the equilibration process, then it would be much simpler 
to describe. The formalism of thermal field theory for QCD would then be 
applicable. If the system is, on the other hand, still out of equilibrium, this 
much more complicated situation would require a phcnomcnological approach 
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as the out-of-equilibrium thermal field theory is not yet fully developed. An- 
other reason is the state of equilibrium of the system has important effects on 
the particle emissions or their production rates. We will show this effect on 
photon production from the plasma. Before we do that, the parton plasma as 
a rather unique many-body system will be shown first. 

2 Equilibration of a parton plasma — fixed vs. evolving a s 

To study the equilibration of the parton plasma, we must have a means of 
performing the time evolution of such a many-body system. To do this, we 
choose some plausible initial conditions from event generator such as HIJ1NG. 
Then we build the time evolution equations from Boltzmann equation 

(l +w *-IO /(p ' T)=cr(p ' T) - (i) 

The collision terms C(p, r) will be constructed in two ways in order to close 
the equations. The first is to use the relaxation time approximation and the 
second is by explicit construction from QCD matrix elements and the particle 
distributions f g , f q and fg. We include all binary interactions at the leading 
order in the renormalized strong coupling 

99 < — ► 99 19 < — > 99 19 * — ► 19 

11 < — ► 11 11 < — > 11 11 * — ► 11 (2) 

and in order to have non-conservation of the number of partons, which one 
would expect from interactions generated by the non-Abelian QCD Lagrangian, 
we include also 

99 < — > 999 99 < — > 11 • (3) 

Combining all these ingredients, we can solve for the distributions, which are 
functions of time and from which most information of the plasma can be ob- 
tained and therefore the time evolution is known. Of course, it is still necessary 
to choose a most appropriate value of the coupling to evaluate the collision 
terms. For this purpose, it is customary to assume an average momentum 
transfer of 2 GeV for the parton interactions in the plasma. This translates 
into an a s = 0.3 for Aqcd = 200 McV. With this choice, the system can be 
evolved in time and one can check for the state of equilibrium. Concerning the 
latter, there are two aspects. One is the balance of the partonic composition in 
the plasma or parton chemical equilibration. The degree of chemical equilib- 
rium is most simply parametrized by the quantity known as the fugacity which 
we denote here by l g and l q for gluons and quarks respectively. The plasma 
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Figure 1: Variation of the gluon and quark fugacities, l g and l q , with time. The solid (long 
dashed) lines are the results of fixed (evolving) a s ■ 



is in chemical equilibrium if / = 1.0 or very far from equilibrium if I = 0.0. 
In Fig. [j], the results of the change of the fugacities with time at LHC and at 
RHIC are shown in solid lines. As can be seen, chemical equilibration for gluon 
is very good at LHC and good at RHIC. On the contrary, that for quarks at 
both colliders is not so good and even poor at RHIC. 

As for thermalization or kinetic equilibration, there is no one parameter 
as in chemical equilibration to parametrize the degree of thermalization, so 
an indicator of the degree of thermalization must be found. Usually, this is 
done by examining the slope of the particle pr-spectra or that of the log of the 
particle distributions. In our case, we choose instead to plot the longitudinal 
to transverse pressure ratios because the pressure must be the same in all 
directions when the system is in kinetic equilibrium. So the closer are these 
ratios to unity, the nearer is the parton system to full thermal equilibrium. 
In Fig. ||, we plot these ratios in solid lines for gluon and quark separately. 
The particular shape of the curves has to do with our rather special initial 
conditions, which we take to be momentarily thermalized at the beginning. 
Therefore these curves all start with a ratio of 1.0. As the expansion sets in, 
the system is driven out from the initial transient kinetically equilibrated state 
and the pressures are no longer isotropic. The system naturally responds to the 
expansion and increases the collision rate. At some point, the latter dominates 
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Figure 2: Variation of the ratios of longitudinal to transverse pressure for gluon and quark 
with time. The solid (long dashed) lines are the results of fixed (evolving) a s . 



over the expansion and the curves all approach unity again after the dip. We 
are more interested in the final approach to kinetic equilibrium. As can be 
seen, the gluons are better equilibrated than the quarks, which are again not 
so good. 

So after we used some reasonable initial conditions, an essentially pertur- 
bative time evolution scheme and an appropriate value of the coupling, the 
state of equilibrium of the parton phase is reasonably good for gluons but not 
so good for the quarks. One wonders whether it is possible to do better than 
this. Evidently, the closer is the plasma to full equilbrium the better. In fact, 
the answer is yes. Improvement is possible because what we have done so 
far, although apparently reasonable, is not entirely correct for a time-evolving 
system. It is this last aspect, which is unusual within the framework of pei±ujj= 
bative calculations of strong interactions, that one could easily overlooktmoQ 
and indeed has been the case until recently a. 

In fact, during the time evolution of the system that produced the previous 
results, the average parton energies dropped significantly by at least 1.0 GeV 
and so the system underwent substantial changes. It is therefore very unlikely 
that the average momentum transfer in an average parton collision can stay 
fixed at around 2.0 GeV throughout. If so, a way must be derived to replace 
the fixed value of a s = 0.3 used to obtain the previous results. A simple recipe 
naturally suggests itself, which is to take the average momentum transfer (Q) 
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to be given by the average parton energy (E). Then one substitutes this into 
the one-loop running coupling expression 

4-7T 

- f>,MAEriK hcD) < 4 » 

to obtain an interaction strength entirely determined by the system. As a 
consequence, this strong coupling will evolve with the plasma and as such this 
new approach is more self-consistent. 

The new results of the time evolution scheme with the self-consistent cou- 
pling are the long dash curves in Fig. ^ and Fig. ^. The dashed curves in Fig. ^ 
and those in Fig. ^ during the final approach to equilibrium all rise faster with 
time than when the plasma is evolved with a fixed a s (solid line). Not only is 
that the case, the final state of equilibrium is markedly better than before at 
both colliders. So a self-consistent time evolution speeds up and improves the 
equilibration of the parton plasma. There are other associated effects on the 
plasma that go along with the aboves but we shall not discuss them here. We 
refer the reader tou. 



3 SU(3) non-Abelian parton plasma is an unique many-body sys- 
tem 

In the previous section, we saw how the evolving coupling improved the equili- 
bration of the parton plasma. This is but a manifestation of an unique property 
of a many-body system governed by the QCD Lagrangian. In the Boltzmann 
equation, it is the collision terms that are responsible for bringing the plasma 
into equilibrium. The collision terms are all made up of the difference between 
the reactions going one way and in the reverse direction. In an ordinary many- 
body system, as equilibrium is approached, the forward and backward reaction 
rates get closer and closer to each other and the equilibration rate will become 
slower and slower as a result. On the contrary in a parton plasma, because the 
interactions are mediated via SU(3) non-Abelian gauge bosons, the interaction 
strength varies with the scale of the interactions. The collision terms of such a 
system are therefore each given by a certain power of a varying a s multiplied 
by the difference of the rate of the forward and backward reaction. As equi- 
librium is approached, the reaction rates in both directions tend to equalize, 
however, the equilibration rate of a parton plasma does not slow down, unlike 
an ordinary plasma, because the increasing strength of the interactions caused 
by the lowering of the average energy of the system is able to more than com- 
pensate for the otherwise equalization of the forward and backward reactions 
or in other words, the slowing down of the equalibration rate. The manifes- 
tation of this effect can be most easily seen on the variation of the collision 
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Figure 3: Photon production at LHC and at RHIC from a parton plasma not in equilibrium. 
The dotted (dashed) lines are from Compton scattering (annihilation) contribution. The 
solid lines on the top are the total contributions. The contribution from Compton scattering 
does not dominate over that from annihilation at LHC. 



time with time (seeEl). Therefore the parton plasma exhibits a most unusual 
feature of the tendency to equilibrate faster and faster. Obviously, this cannot 
continue forever and at some point the equilibration of the quarks and gluons 
will be interrupted by the start of the deconfinement phase transition. 



4 Out-of-equilibrium effect on photon production 

We have shown the time evolution of a parton plasma at the two colliders. 
All the information of the system is contained in the particle distributions 
and from which we learned about the plasma's state of equilibrium. With 
the knowledge of the momentum distributions of the particles, one can work 
out the rates of particle production such as photons, dileptons etc.. Also our 
method of evolving the plasma in time, imposed no equilibrium requirement on 
any aspects of equilibration and hence the form of the distributions. As such, a 
direct comparison of particle productions of a plasma that is out of equilbrium 
with an equilibrated one can be done. We will show an out-of-equilibrium 
effect on photon production and explain its importance. 

Photon production from the plasma is dominantly through Compton scat- 
tering and annihilation contribution. The total production rate from these two 
contributions is given by 

P d 7 N 1 f d 3 kx d 3 k 2 d 3 k 3 4 (4) 

x |2/ 3 (fc 1 ,r)/ 9 (fc 2 ,r)(l-/ 9 (fc 3 ,r))|X 9 „ 7 | 2 

+f q (k u r)f l} (k 2 ,T)(l + f g (k 3 ,T))\M qq ^ g - l \ 2 ) . (5) 
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In ref. B, a kinetically equalibrated plasma was studied and the photon 
rate calculated. It was shown in their Fig. 8 and 9 that the contribution from 
Compton scattering was far more important than the annihilation contribution 
both at LHC and at RHIC energies. If one now compares these figures with 
the one we have here in Fig. one can immediately see that in our plots, 
Compton scattering contribution does not dominate over that of annihilation. 
In fact at higher px at LHC, annihilation contribution is the more important. 
At RHIC in the same px range, although Compton scattering is larger, it is not 
by the same large amount as shown in ref. El. The reason for these differences 
between our photon production and thosejn ref. □ is because our plasma is 
not in kinetic equilibrium but that in ref. a was. We will explain this in the 
next paragraph. 

In Eq. ([5]), the Compton and annihilation contribution to photon yield 
differ by various factors such as the interaction matrix elements, the Pauli 
blocking and Bose-Einstein stimulated emission, and one of the initial distri- 
butions. In fact, the matrix elements do not differ significantly numerically 
and the Pauli blocking and Bose-Einstein stimulated emission also provide 
no great suppression or enhancement. The source responsible for the differ- 
ence seen in the two different sets of figures comes actually from the gluon 
distribution in Compton scattering and the initial quark distribution in the 
annihilation contribution. To explain this more clearly, we use the observation 
that the difference in the figures is most prominent for larger px photons — 
annihilation contribution is actually larger than Compton scattering at LHC 
in Fig. H — so we can concentrate in this pt range. Emission of larger pt 
photons requires more energetic incoming partons, this let us simplify the ex- 
planation by assuming Boltzmann form for the incoming particle distributions 
/ ~ I exp{— k° /T}, while for the outgoing ones, we will replace them by unity. 
In reality, Pauli blocking and Bose-Einstein stimulated emission do have some 
effects but they are not so important in comparison with the one that we are 
going to describe so it is best to leave them out. Since we are comparing the 
two contributions, we form the ratio 

after taking \M. gq tql \ ~ \M-qq ^l- Because of the initial gluon dominance 

and the stronger interaction amongst gluons due to colour, l q /l g << 1 always 
and therefore if the exponential factor is not present, Compton scattering in- 
deed dominates over annihilation contribution. In a kinetically equilibrated 
parton plasma, the system has only one temperature T g — T q , so the expo- 
nential factor is unity and we get the results shown in Fig. 8 and 9 of ref. 
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In our present case, the plasma is not in kinetic equilibrium but we can 
most simply consider them as a mixed fluid of gluons and quark- antiquarks 
at different temperatures. The gluon temperature is higher initially but will 
eventually be lower than that of the quarks because of gluon multiplication and 
their conversion into quark-antiquark pairs. When that happens, for a photon 
with a high enough px, the incoming parton energies will be sufficiently large 
that the product in Eq. (||) will be larger than unity and the photon rate from 
annihilation will be larger than that from Compton scattering. This is what 
one sees in our result in Fig. ||| at LHC. In the same figure at RHIC, there is 
not sufficient time for the two temperatures to be in the correct range from 
each other so that the exponential factor is able to compensate for the fugacity 
ratio. The result is Compton scattering remains larger than the annihilation 
contribution iut it is definitely by a lesser amount than that shown in Fig. 8 
and 9 of ref. a. 

We see that the out-of-equilibrium effect tends to enhance the annihilation 
contribution at higher pt and hence the total photon yield in that pt range. 
Therefore photon production of partonic origin should have a better chance to 
compete with the hard photons from the initial collisions and those fragmented 
off minijets. The window for observing photons from deconfined matter will 
be widened as a result. 
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